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Abstract 

We have solved Cosmological Gravitational Wave(GW) equation in the frame work of Gen- 
eralised Brans-Dicke(GBD) theory for all epochs of the Universe. The solutions are expresed in 
terms of the present value of the Brans-Dicke coupling parameter u ((/)). It is seen that the solu- 
tions represent travelling growing modes for negative values of loq for all epochs of the Universe. 
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1 INTRODUCTION 

General Relativity (GR) predicts that perturbations of a background spacetime generate grav- 
ity waves. The general results concerning gravitational waves have been investigated by numerous 
authors long ago[l].In the context of GR cosmologies, Grishchuk [2] has shown that the varying 
gravitational field of the expanding Universe would amplify zero -point fluctuations, and lead to 
the formation of relic gravitons .Gravitons interact very weakly with matter and other fields, and 
only at energies of order of the planck enegy( 10 19 GeV) will the gravitational interaction become 
as strong as the electromagnetetic interaction. Therefore, the relic gravitons must have decoupled 
from the dynamics of the Universe at extremely early times, and offer a unique way to probe the 
structure and evolution of the early Universe. Even without detecting them, we can still get use- 
ful information by studying the consequences of their presence in several processes of the cosmic 
history ;for example, primordial nucleosynthesis and the anisotropy they would induce in the mi- 
crowave background .Thus, the idea of using gravitational radiation as probe of the early Universe 
has turned into a realistic possibility. The subject of gravitational radiation within Brans-Dicke 
(BD) theory was carried out by Wagoner [3] This was concerned with gravitational waves generated 
in the weak field limit. Barrow[4]has derived the cosmological gravitational wave equation in some 
more general scalar tensor theories of gravity and solved the wave equation for universe in BD 
theory where the coupling function is a constant. However, in the generalised Brans-Dicke(GBD) 
theory of gravity where the coupling function uj = uj ((f)), the gravitational wave equation is solved 
by Barrow for only vacuum and radiation dominated epochs of the Universe for a particular choice 
of the coupling function a; ((/>). Thus the work seems to be incomplete in GBD theory of gravity.Here 
we investigate the cosmological solutions of the gravitational wave equation for all epochs of the 
Universe in GBD theory by allowing the coupling parameter to run with scalar field <fi following the 
approach adopted in our recent work [5]. The solutions are expresed in terms of present value of the 
coupling parameter u)((j)) i.e,in terms of ojq. An outline of this paper is as follows. In the section II 
we have derived the solutions. The solutions are expresed in terms of lvq in section III for all epochs 
of the Universe and shown that for negative values of ujq solutions represent growing modes. .We 
conclude with section IV. 



2 WAVE EQUATION 



The gravitational field equations are given below [6] 
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When perturbation is applied to above equations,they become[4,7] 
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Where the perturbed metric is 

= + V" where,%„ = h^ v 

leading to the Ricci tensor, Rqq = -^oo + SRoo- 

The energy-momentum tensor is given by , T 00 = T 00 + (5T 00 = T 00 + Sp, 
The trace of energy-momentum tensor is ,T = T + 5T = T + 5p — 35P. 

For the d'Alembertian of the scalar field we find , 5u<j) = 54> + aah ij <j) - ^h^ + 3^5<p - %r5<f>. 
We adopt the synchronous gauge, ho^ = 
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The time time component of perturbed Ricci tensor is[8], 5Rqo 
When <^>=constant,all the perturbed equations will reduce to perturbed equation of Eienstein. The 
perturbations to the metric that represent weak gravitational waves can be expressed as [9] 
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Where the functions Y k (t) and Qj(k, x)satisfy 
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aplacian operator, 
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In the above equations, a(t) is the scale factor, i,j = 1,2,3; V 2 is the spatial 
an over dot indicates derivatives with respect to the synchronous cosmic time ,t, and k = \k\ = 
with k being the comoving wave vector. The eigen functions Qj(k,x) can be written in terms of 
plane -wave solutions exp(±ik.x) times a constant polarization tensor and using the field equations 
in the zero curvature FRW model we can simplify (3) and recast it as [4] 
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This is the Gravitational wave equation in more general scalar-tensor theory of gravity. 



3 SOLTUIONS IN GBD THEORY 



The generalised Brans-Dicke scalar-tensor theory of gravity is characterised by the generalisation 
of the coupling function w to a function of eft; and U = ^ = O.We introduce the conformal time ij 
defined by dt = adrj 

Each component hff of the gravitational wave perturbations can be written as[ 4] 
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Here R(r]) = a(r/)0 1 / 2 (?7) and 
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with the prime indicating derivatives with respect to ij. 
When Lo((j))= constant, eqn.(6) becomes, fj,"(k,rj) + 
which represents Gravitational wave equation in Brans-Dicke scalar-tensor theory of gravity.The 
solution given by eqn(6) is formally identical to the corresponding one obtained in the general rela- 
tivistic case, but with the quantity R(r]) defined by eqn(5) replacing the scale factor a(rj). Equation 
(6) describes an oscillator with varying frequency or a Schrodinger equation with potential 

_R» dm/d<P <j)> 2 
R ~ l ~ 2a;+3 <f> 

for a particle with energy k 2 and with the variable rj playing the role of a spatial co-ordinate. We 
solve eqn(6) for different epochs of Universe,with equation of state P = jp, — 1 < 7 < 1. 
The scale factor and the scalar field evolve as power laws of the cosmic time 
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For a / l,the translation of eqn(7) and (8) into the conformal time yields 
a(v) = C {r]-rg) q 
0(v) =<ro(v- VoY 
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The wave equation (6)takes the form 
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We take 2to + 3 = (2lo + 3)<f> n 
This leads to 2^ = (2lu + 3)n(f> 
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This reduces equation (9)to 
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Hence,we get 
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For very large wave numbers or,more precisely, for waves such that 
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fc 2 > ^(r/), where F(r?) = |J2"/i2 + ^^g^- 
eqn(10) will have purely oscillatory solutions, and eqn(5) takes the form 



where Ci and C 2 are arbitrary constants and the third constant 770 is same as before.We see that 
the effect of the scalar field <fi is for the scale of amplitude decrease to be set by the function R(rj) 
instead of by a(rj). 



In the opposite regime of very small wave numbers, where k 2 <C V(rj), 
we obtain, hf (rj, x) = [d + C 2 J iR- 2 {r/)dr/] Qj(k, x) 

Note , however, that ,for bell-shaped potentials, the amplification coefficient will scale as R(r) f) / R{rfi) 
instead of a(rj f) / a(i]i) . 

The general solution of eqn(10) can be written as [10] 
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Here J m and Y m are Bessel and Neumann's functions respectively ,X(rj) = k(rj — fjo) and 
m = s' - 1/2 
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We see that in order to represent a physical mode the power m has to be real and greater than zero. 
The difference between the solution (11) and the corresponding one in GR lies in the presence of 
the coefficient (3 in the parameters s' and m.Setting (3 = we recover the GR solution as expected, 
since this corresponds to the requirement that 4> ,and hence G, remains constant. The general 
cosmological solutions in the framework of GBD theory when matter satisfies the equation of state 



P = 7/9 with — 1 < 7 < 1 a constant,were derived by Sahoo and Singh [11]. An important feature of 
these solutions is that they are expressed in terms of the present values of the coupling parameter 
a;, for vacuum,radiation and matter dominated epochs of the Universe.We write gravitational wave 
perturbations in different epochs of the Universe by cosidering these solutions in terms of time t. 



3.1 Vacuum epoch 

For this epoch a,/3,n are [11] 

a = ZJ^2^ = n = "(2^0 + 5 ) 



1 +"0 



As,dr] = f ,we find,!? « (f±^)i 2 +-o 
And we take,(r? - rjo) « (j^)t^ 



Hence, 



1 +"0 



r=0 



+ 



It is clear from the above expression that,o>o 7^ —2 which is in agreement with the conclusions 
obtained in our earlier investigation [11]. In order to represent \i a growing modes,the power of t 
should be real and greater than zero. This gives following conditions 

(i) wo > -1 

(ii) (5 + l)(m + 2r) > -\ 

(iii) m should be real. 

In this case, the order of the Bessel function is m = iiwl + 2 ^ Q+9 



As m should be real,we get that ujo > — | which is automatically satisfied if the condition (i) 
i.e,a>o > — 1 is valid. Thus,we conclude that for gravitational waves to grow ujq should take values 
greater than -1. 



3.2 Radiation epoch 



For this epoch ce,/3 and n are [11] 
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Thus, for gravitational waves to grow during the radiation epoch the above solution should satisfy 



the following conditions, 

(i) uj > -3 

(ii) (q+ l)(m + 2r) > -\ 

(iii) m should be real. 

Using the values of a,f3 and n in the expression for m ,we get m = - As m should 

be real we get ujq > — | which automatically satisfies constraint(i).Thus,for gravitational waves to 
grow wo should take values greater than — |. 

^ = i MM) C«(M) 

hij = fd 3 k hfj(k,t) 

This represents the perturbation in Gravitational field during the radiation epoch. 
3.3 Matter epoch 



During this epoch a, (3, and n are [11] 
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It is clear from the above expression that,a>o 7^ —2 .Again this result agrees with our earlier conclu- 
sions^!]. For growing modes,the power of t should be real and greater than zero. This gives following 
conditions 
(iVo > -2, 

(ii)(g + l)(m + 2r) > -\ 
(hi) m should be real. 

In this case,the order of the Bessel function,m,is m = ± \ yT + (2+w$ • 
As m should be real,we get ooq > —3. 

Thus,for perturbation to grow u)q should take values greater than -2. Finally, we get = ^ n{k,t) Cij (k, t) 
and 

h ij = fcPk hlj{k,t) 

This represents the perturbation in gravitational field during the matter epochs. 

4 CONCLUSIONS 

We have solved the gravitational wave equation for a general class of scalar tensor theories of 
gravity. The solutions are expressed in terms of the present value of the Brans-Dicke coupling pa- 
rameter i.e,u;o.The solutions represent growing modes when the parameter ujq takes negative values 
particularly greater than -l.This negative values of the parameter are required for structure for- 
mation, cosmic accleration,superluminal expansion and at the same time a significant amplification 



of perturbations [7, 12] .It is also seen that,unlike in GR, there can be amplification of gravitational 
waves in a radiation -dominated universe. These new features are consequence of the more general 
scalar tensor theories of gravity which in special cases reduces to Brans-Dicke thory of gravity 
and GR. In a subsequent paper we will use the results derived here to obtain the spectrum of the 
relic gravitons in this model.This will enble us to constrain some of the parameters appearing in 
the scalar -tensor theories. In particular, it could impose a further restriction on the value of the 
coupling parameter at early times. A small tu(<f)) is required to have satisfactory nucleation process 
in some inflationary cosmologies. 
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